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Abstract 

The  radiative  losses  due  to  free-free  transitions,  free-bound 

transitions,  and  permitted  and  forbidden  line  emission  have  been  evaluated 

for  all  ions  which  are  present  in  significant  amounts  in  the  extended  solar 

corona  from  4R  to  215R  .  Most  of  the  radiative  losses  occur  interior  to 
e  a 

15R  ,  and  they  are  dominated  by  line  emission  in  the  extreme  ultraviolet 
and  x-ray  wavelengths.  In  general,  the  radiative  losses  are  found  to  be 
unimportant,  and  the  consequences  of  this  result  are  discussed  in  terms  of 
the  semiempirical  model  recently  presented  by  Brandt  and  Mi  chie  (1962). 
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1:  Introduction 


The  determination  of  the  radiative  losses  in  the  solar  corona  and 
extended  solar  atmosphere  is  of  importance  in  discussions  of  the  energy 
balance  and  structure  of  the  corona.  Earlier  estimates  of  these  losses  have 
been  made  (Woolley  and  Allen  1948;  Elwert  1954;  Zirin  1957;  Allen  1961;  Lvist 
and  Zirin  I960;  De  Jager  and  Kuperus  1961).  However  the  observational 
values  of  many  of  the  physical  quantities  in  the  corona  and  interplanetary 
medium  have  not  been  made  certain  as  yet.  As  work  progresses  and  more 
observational  material  is  obtained,  the  calculation  of  the  radiative  losses 
which  are  based  on  this  data  must  continually  be  revised. 

In  a  recent  paper,  Brandt  and  Michie  (1962,  Paper  VII)  obtain  a 
theoretical  expression  for  an  expanding  solar  corona  which  they  find  to  be 
reasonably  consistent  with  empirical  evidence  (Brandt  1962a);  recently, 
the  values  of  the  expansion  velocity  and  density  found  by  Brandt  and  Michie 
(1962)  have  been  confirmed  by  the  results  obtained  from  Mariner-II  by 
Neugebauer  and  Snyder  (1962).  In  such  an  expanding  corona,  and 

are,  respectively,  the  total  energy  per  second  and  the  total  mass 
per  second  which  cross  a  spherical  surface  at  a  distance  r  from  the  sun, 
then  at  a  point  beyond  which  energy  deposition  and  radiative  losses  can  be 
neglected,  the  ratio  will  be  a  constant  independent  of  r.  Brandt 

and  Michie  have  evaluated  this  constant  at  4  solar  radii,  thereby  introducing 
the  assumption  that  radiative  losses  and  energy  deposition  are  negligible 
beyond  that  point.  It  is  the  purpose  of  this  paper  to  test  the  validity  of  the 
assumption  with  regard  to  the  radiative  losses. 

2.  Preliminary  Results 

In  order  to  find  the  rate  of  energy  lost  by  radiation  in  the  entire 
region  between  4  solar  radii  and  the  orbit  of  the  earth,  we  must  evaluate 
an  integral  of  the  form 


-3- 


T  = 


21 5R 


€(r)r2dr 


erg/ sec 


(1) 


where  «  (r)  is  the  rate  of  energy  loss  per  cm  by  all  radiative  processes. 

The  quantity  €(r)  may  be  broken  up  into  three  contributions;  continuous 
emission  due  to  free-free  transitions,  continuous  emission  due  to  free- 
bound  transitions,  and  line  emission.  These  contributions  must  be 
evaluated  for  all  ions  which  are  present  in  significant  amounts  in  the' 
extended  solar  corona. 

The  abundance  of  the  various  ions  is  determined  first  from  an 
assumption  concerning  the  abundances  of  the  elements  in  the  region  vinder 
consideration  and  second  from  a  calculation  of  the  state  of  ionization  of  the 
various  elements.  We  base  our  calciolations  upon  the  relative  abundances  of 
the  elements  as  given  by  Goldbert,  Muller,  and  Alley  (I960);  the  elements 
considered  and  their  relative  abundances  are  indicated  in  Table  1. 

To  determine  the  state  of  ionization  of  the  various  elements  we 
note  that  the  ionization  equilibrium  in  the  corona  is  established  by  a  balance 
between  radiative  recombinations  and  excitation  by  collisions  with  electrons, 
under  the  assumption  that  photo- ionizations  and  three-body  recombinations 
can  be  neglected  (Woolley  and  Allen  1948).  However  the  results  of  Brandt 
and  Michie  (1962)  indicate  that  the  gas  of  the  corona  is  expanding  outwards 
with  velocities  ranging  from  about  8  km/sec  at  4  solar  radii  to  300  km/sec 
at  the  orbit  of  the  earth.  Because  of  this  expansion,  the  material  at  a 
given  point  may  not  have  had  time  to  reach  secular  equilibritim  between  the 
various  stages  of  ionization. 

Consider  two  points,  A  and  B.  The  time  for  material  to  travel 
from  A  to  B  is 


t. 

J 


~  ^ab'^'^^^ab 


(2) 
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where  is  the  distance  from  A  to  B  and  <  v  >  is  the  average  velocity 

between  A  and  B.  If  the  ions  in  stage  of  ionization  p  recombine  and  no  new 
ions  in  stage  p  are  produced,  then  the  time  for  the  number  of  ions  in  stage  p 
to  drop  to  1/e  of  its  original  value  is 


t 

r 


=  <  a 


Ne> 


-1 


(3) 


where  is  the  electron  density  and  Q  is  the  recombination  coefficient  to 
stage  p-1.  The  qixantity  t^  then  is  a  characteristic  time  for  recombination 
of  ions  when  the  material  moves  to  a  region  favoring  lower  ionization.  A 
sufficient  condition  that  the  relative  number  of  ions  in  a  particular  stage 
not  change  when  the  material  moves  from  A  to  B  is 


t. 

J 


<< 


t 

r 


(4) 


We  find,  using  the  estimates  of  expansion  velocity  by  Brandt  and  Michie, 

that  for  A  =  2R^  and  B  =  4R^,  t.  =  t^/8;  and  for  A  =  4Rg  and  B  =  ISRg, 

t.  =  t  /20  for  an  ion  of  charge  Z  =  14;  furthermore  t.  will  be  even  smaller 
J  r  »  j 

compared  to  t^  for  ions  of  a  smaller  charge.  Therefore  we  conclude  that 
the  state  of  ionization  of  the  elements  in  the  outer  corona  beyond  4  solar 
radii  is  the  same  as  it  is  in  the  corona  within  4  solar  radii. 

In  the  corona  we  may  assume  that  secular  equilibrium  holds, 
that  is,  at  a  given  point  the  rate  of  formation  of  ions  of  stage  p  equals  the 
rate  of  disappearance  of  ions  from  stage  p.  An  ion  in  stage  p  can  be 
created  through  recombination  from  stage  p  +  1  or  through  ionization  from 
stage  p-1.  Under  the  assumption  that  three-body  recombinations  and 
photoionizations  can  be  neglected,  we  arrive  at  the  secular  equilibrium 
equation 

^e^p^®p-»p+l  “p^p-1^  “  ^^p-l®p-l-4  p  ^p+l®p+l— »  p^  ^e 


(5) 
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where  N  represents  the  electron  density,  N  , 
e  p 


Np  j,  and  represents 


the  ion  densities  of  stages  p,  p-1,  p+1,  respectively,  S.  .  represents 
the  collisional  ionization  coefficient  for  stage  i,  and  represents  the 

recombination  coefficient  onto  stage  j.  If  stage  p  corresponds  to  a 
completely  ionized  atom,  equation  (5)  reduces  to 


Na.,=N,S, 

p  p-»p-l  p-1  p-l-»p  (6) 

Now,  by  replacing  p  by  p-1  in  equation  (5)  and  by  substituting  equation  (6) 
into  the  result,  we  obtain  another  equation  of  the  same  form  as  (6).  Thus 
we  find  that  for  each  stage  of  ionization  the  following  relation  (due  to  Woolley 
and  Allen  1948)  holds: 


JL  =  fp-tP+L 

“p+l-#p 


where  (Woolley  and  Allen  1948,  Allen  1955b) 


S  =  3  X  10*® 

p-^p+1  1 


-llbOOXi 


exp 


(8) 


Q  =  1.  5  X  10'^^ 


(9) 


Z  is  the  ionic  charge  of  stage  p+1,  is  the  ionization  potential  of  stage  p 

in  ev,  and  x  is  the  number  of  ions  in  the  p+1  stage  divided  by  the  total 
number  of  ions  in  the  p  and  p+1  stages. 

From  equation  (7)  we  calculate  the  ion  densities  of  the  more 
abundant  species  relative  to  hydrogen,  assuming  a  temperature  of  10®*^K 
and  using  the  cosmic  abundances  from  Table  1.  The  results  are  shown  in 
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Table  2.  The  method  ia  not  used  for  the  elements  hydrogen  and  helitim;  we 
assume  that  both  elements  are  essentially  completely  ionized  in  the  region 
imder  consideration. 

The  calculation  of  radiative  loss  also  requires  a  knowledge  of  the 
temperature  and  electron  density  distribution  in  the  outer  corona.  The 
temperature  distribution  near  the  sun  can  be  obtained  from  the  observed 
densities  by  use  of  the  assumption  of  hydrostatic  equilibrium  (Pottasch 
I960,  Chapman  1961).  We  use  two  different  temperature  distributions, 
both  of  which  agree  with  the  Chapman  values  between  4  and  5  solar  radii. 
Distribution  (A)  uses  a  probable  upper  limit  of  100,  000°K  for  the  tempera¬ 
ture  at  1  A.  U.  ,  while  distribution  (B)  uses  a  probable  lower  limit  of 
30,  000°K  at  that  point  (Chamberlain  1961).  The  intermediate  values  are 
obtained  by  interpolation. 

Two  distributions  of  electron  density  are  obtained  in  a  similar 

manner.  Between  4  and  6  solar  radii  the  values  in  both  cases  are  those 

observed  by  Blackwell  (1956,  see  also  van  de  Hulst,  1950).  Near  the  orbit 

of  the  earth  the  value  =  1  (Brandt  1962a,  Ingham  1961)  is  used  in 

distribution  (a)  while  in  distribution  (b)  the  value  N  =  4  is  used  as  a 

e 

reasonable  upper  limit.  The  intermediate  values  are  again  interpolated. 
Since  we  take  the  ratio  of  abundances  of  hydrogen  to  helium  as  10:1  and  since 
both  are  completely  ionized,  the  density  of  hydrogen  ions  is  =  0.  83N^. 

In  our  calculations  we  shall  assume  a  Maxwellian  velocity 
distribution,  which  is  a  reasonable  approximation  at  points  close  to  the  sun 
but  which  breaks  down  at  large  distances.  Fortunately  the  calculations  show 
that  most  of  the  radiative  loss  occurs  at  distances  of  less  than  15  solar  radii, 
therefore  the  introduction  of  this  assumption  should  not  result  in  appreciable 
error. 


3.  The  Method  of  Calculation 

We  now  turn  to  the  actual  calculation  of  the  quantityc  (r).  The 
intensity  of  the  free-free  emissions  is  given  by  (Allen  1955a) 
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e  ,,  =  1.  44  X  10”^^  N  N.  ergs/cm^sec 

11  e  1  " 


(10) 


where  g  is  of  the  order  unity,  Z  is  the  ionic  charge,  N  and  N.  are  the 

36  1 

numbers  of  electrons  and  ions,  respectively,  per  cm  .  Using  the  faict 
that  for  a  given  ion 


N.  =  0.  83  N  p. 
1  e^i 


(11) 


where  p.  is  the  abundance  of  the  ion  relative  to  hydrogen,  we  obtain 

€  ,,  =  1. 19  X  10'^^  ^  ^  Z.^p.  ergs/cm^sec  (12) 

11  e  .  1  1 

1 


where  the  summation  is  over  all  the  ions  of  interest.  Calculation  of  the 

E2 

Z.  p.  gives  a  result  of  1.  4  for  hydrogen  and  helium  and  0.  076 
for  all  the  other  elements.  Thus  the  final  expressions  become 


€  „(H  +  He)  =  1.  67  X  lO'^^  N  ^ 
ff  e 

t  ^^(total)  =  1.  76  X  10'^^ 


ergs - 

cm  sec 


(13) 


The  rate  of  free-bound  emission  is  derived  by  Menzel  (1937);  the 
total  energy  emitted  in  transitions  from  the  continuum  to  a  level  n  of  a 
hydrogen-like  atom  is 


^  =  N^N^kKZ^gbj^/n^T^^^  ergs/cm^sec 


€ 


(14) 
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where  bj^  indicates  tlie  departure  from  a  Mascwrellian  velocity  distribution. 
Summing  over  n,  evaluating  the  constants,  setting  b^^  =  1,  and  making  use 
of  equation  (11)  we  have 


€  ^  =  4.  52  X  10~^^  ^  ergs/cm^sec  (15) 

This  formula  is  applicable  to  the  calculation  of  the  free-bound  energy 
loss  for  H  II  and  He  III  only;  for  these  we  find  that  ~ 

the  heavy  elements  we  use  equation  (21)  below. 

In  the  calculation  of  the  bound-bound  energy  loss  we  first  note 
that  at  the  dilute  radiation  densities  and  low  particle  densities  vinder 
consideration,  the  times  for  particle  or  photon  excitation  or  de- excitation 
of  bound  electrons  are  much  greater  than  spontaneous  radiation  times 
(<vlO“  sec).  Thus  an  electron  in  an  excited  state  will,  in  general,  cascade 
to  the  ground  state  before  any  other  process  can  occur.  Therefore  part  of 
the  bound-bound  emission  will  be  due  to  electrons  which  cascade  to  the 
ground  state  after  they  arrive  at  levels  n  =  2  or  higher  by  processes  of  ^ 
recombination.  To  find  the  energy  from  this  source  we  multiply  the  total 
number  of  recombinations  to  each  level  by  the  energy  difference  between 
that  level  and  the  ground  level.  The  recombination  coefficient  to  the  nth 
level  of  a  hydrogen-like  atom  is  given  by  Cillie  (1932)  as 


a  =  3.  26  X  10'^n‘^T‘^'^^Z^  exp 
n 


^n/kT/  E^(J::n/kT)  (16) 


where 


/•fc  e”*" 

Ej(x)  =  /  — — dt 


(17) 


and  where  X  ^  for  a  hydrogen-like  atom  is  given  by 
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157900Z^ 

XnjV.T  = - X -  =  .  (18) 

n  T 


The  energy  is  then 


bb(l) 


“n^i^e  ergs/cm^sec 


(19) 


where  ^  is  the  ionization  potential.  This  formula  is  used  in  the 
calculations  in  the  case  of  hydrogen  and  helium.  For  the  heavy  elements 
we  estimate  the  energy  loss  due  to  free-bound  transitions  and  subsequent 
cascades  to  the  ground  level  by  using  the  generalized  recombination 
coefficient  of  a  non-hydrogen-like  atom,  given  by  Allen  (1955b): 

Q  =  1.  5  X  lO'^^Z^T"^''^  .  (20) 


The  energy  loss  becomes 

*  fb  +  bb(l )  17  elements^  =a  N^N.(r  +  2 


(21) 


where  (  is  the  average  kinetic  energy  of  the  recombining  electron,  assumed 
to  be  100  ev.  This  formula  is  used  for  all  ions  listed  in  Table  2, 

The  second  source  of  bound  -bound  emission  is  collisional  excitation 
by  electrons  from  the  ground  level  followed  by  cascade  back  to  the  groxmd 
level.  For  the  case  of  hydrogen  and  helium  we  follow  the  method  of  Zirin 
(1957).  To  determine  the  energy  due  to  bound-bound  transitions  he  assximes 
that  ionization- re  combination  equilibrium  holds: 


N(H  I)  X  ionizations=  N(H  II)  x  recombinations. 


(22) 


-1^ 


However  our  calculation  is  based  upon  a  model  of  an  expanding  corona; 
therefore  we  would  not  expect  such  an  equilibrium  to  hold.  Since  hydrogen 
is  completely  ionized  in  the  inner  corona,  the  number  of  H  II  ions  at  any 
point  in  the  outer  corona  is  actually  greater  than  the  number  one  would  expect 
under  conditions  of  equilibrium.  CXir  calculation  of  the  bound- bound  energy 
from  this  source  therefore  represents  an  upper  limit  since  it  is  based  on 
the  equilibrium  assumption. 

The  hydrogen  line  emission  is  given  by  Zirin  (1957,  eq.  4A)i 


recombinations 

f  X  .  .  ,  ■  X  excitations  xX  . 

^  '  ionizations 


(23) 


A  similar  equation  holds  for  helium.  However,  Bethe  {1933a)  shows  that 
for  hydrogen-like  ions,  the  rate  of  collisional  excitation  to  all  levels  and 
the  rate  of  collisional  ionization  are  in  the  ratio  of  about  54:36  for  electron 
energies  of  1000  electron  volts,  which  corresponds  to  an  electron  tempera¬ 
ture  of  10^  °K.  Furthermore,  extrapolation  of  Bethe's  table  indicates  that 
the  ratio  tends  towards  1  for  lower  electron  energies.  Therefore  we  can 
take  as  valid  Zirin's  assumption  that  the  rates  of  collisional  ionization  and 
collisional  excitation  are  roughly  equal  in  the  case  of  hydrogen.  Bethe 
also  shows  (1933b)  that  the  rates  of  collisional  ionization  and  excitation 
both  depend  on  Z  ,  so  the  ratio  is  independent  of  Z  and  the  same  assumption 
can  be  applied  in  the  case  of  helium.  We  use  the  following  expression  in 
the  determination  of  the  bound-bound  energy  loss  from  hydrogen  and  helium, 
based  upon  a  formula  from  Cillie  (1932): 

=  3. 26  X  lO'^N  N.  ^ 

bb(2)  e  1  1  - 3 - 

n 

where  X  ^  is  the  ionization  energy  and  Q  is  given  by  equation  ( 18) . 

We  now  turn  to  the  calculation  of  the  radiative  loss  due  to 
permitted  lines  of  the  heavy  elements.  Because  1/A  »  tj  for  distances 
over  which  T  and  change,  secular  equilibrium  for  excitation  is  valid 


(24) 
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at  any  point,  and  the  number  of  bound  electrons  collisionally  excited  in 

a  given  time  interval  equals  the  number  of  photons  produced.  Following 

3 

the  development  of  Woolley  and  Allen  (1948),  the  energy  loss  per  cm  per 
sec  in  a  given  multiple!  whose  lower  term  is  the  ground  term  is 


*  permitted 


1.  88  X  10"  ^  ^ 

- :pj2 — ^ 


f  p. 
em^i 


[ 


exp(-y)-yEj(y)j 


(25) 


where  y  =  %'/kT,  X'  is  the  excitation  potential  of  the  multiple!,  and 

f  is  the  absolute  value  of  the  emission  f-value  for  the  transition.  The 
em 

f-values  used  in  the  calcvilations  are  those  recently  computed  by  Varsavsky 

(1961). 

A  number  of  important  permitted  lines  fall  in  the  X-ray  region, 
primarily  below  100  A;  f-values  for  these  transitions  are  not  given  by 
Varsavsky.  In  order  to  take  these  into  account  we  make  an  estimate 
based  on  the  calculations  of  Elwert  (1954,  1958).  He  gives  (1958,  eq.  12) 

3 

the  energy  loss  resulting  from  collisional  excitation  per  cm  per  sec  as 


J 


L 


€(X-ray)  =  f^ 


4  X  10 


-19 


(26) 


where  the  summation  is  taken  over  all  the  elements  and  stages  of 
ionization  under  consideration;  x-ray  lines  considered  are  listed  in  Table  2i 
Here,  f^  is  an  uncertainty  factor  of  the  order  unity  and  the  quantities 
include  abundance  factors  as  well  as  collisional  cross  sections.  The  ions, 
the  wavelengths  of  the  lines  considered,  and  the  as  a  function  of 
temperature  are  listed  in  Elwert  (1954),  Table  5.  The  calculation  was  made 
with  our  electron  density  distribution  and  an  assumed  T  =  6  x  10^  °K,  a 
simplification  which  can  introduce  only  negligible  error  since  the  quantity 
turns  out  to  be  only  weakly  dependent  on  temperature.  We  note  that  two 
lines  in  Elwert' s  table,A;40  of  C  V  andA33  of  C  VI,  are  also  included  in  our 
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calculation  of  equation  (25).  A  comparison  of  the  radiative  loss  due  to  these 
lines,  based  on  the  two  independent  methods  of  calculation,  shows  reasonable 
agreement;  therefore  we  feel  justified  in  using  EHwert's  formula  for 
calculating  the  X-ray  emission. 

In  turning  to  the  calculation  of  the  radiative  loss  from  forbidden 
lines,  we  first  note  that  only  photons  arising  from  collisional  excitation 
will  affect  the  flux  balance.  Unfortunately,  the  excitation  cross  section 
has  been  calculated  for  only  one  line,  the  X5303  line  of  Fe  XIV,  in  which 
case  the  excitation  cross  section  is  given  by 

ffiv)  =  0.  78  v‘^  cm^  (27) 


It  follows  that  the  collisional  excitation  rate  from  level  a  to  level  b  is 
(Firor  and  Zirin,  1962) 

C  ,  =  1.  60  X  10"^  (28) 

ab  e 

If  we  assume,  as  do  Firor  and  Zirin,  that  the  excitation  cross  sections 
for  all  the  forbidden  lines  under  consideration  are  about  the  same,  we  may 
estimate  the  emission  due  to  forbidden  lines  by  the  lormula 


*  forbidden  =  1-  ^3  x  10'^  T  p.bv.  ergs/cm^ec  (29) 

1 

In  the  calculations,  all  of  the  forbidden  lines  listed  in  Table  II  of  Woolley 
and  Allen  (1948)  were  considered. 

To  summarize,  the  formulae  used  in  the  calculations  are  (13), 

(15),  (19),  (21),  (24),  (25),  (26),  and  (29).  Substituting  N.  =  0.  83  p.N^ 
into  equations  (19),  (21),  and  (24),  evaluating  ^h  p.  in  equation  (29) 
and  ^  Z  p.(r  +X.)  in  equation  (21),  using^  Y  .  =  27  x  10  ”  in  equation  (26]t 

1  1  Z1 

including  an  estimate  of  the  bound-bound  contribution  from  upper  levels  of 

the  heavy  atoms  in  equation  (25)  (see  Results),  and  combining  terms  where 

possible,  we  arrive  at  the  following  expression  for  the  total  emission  per 
3 

cm  per  sec: 


-13. 


tot 


N  ^UX 


=  1.  76  X  10  'T 


-27^1/2 


["l.  28  X 


10“^®  +  2. 82  X 10“^^  ^ 

heavy 

elements 


+  2.  70  X  f  ^ 

n=l 

i  =H,  He 


exp(Q)Ej(0 ) 

- 3 - 

n 


40 


+  Y  exp(9)Ej(^) 

n=2 

i=H,  He 


n 


where 


y  = 


-jTjr-  and  “X  '  is  the  excitation  potential  for  a  heavy  element 
1579007^ 


=  ionization  potential 

=  excitation  potential  for  H  and  He  II 

p^  =  abundance  of  ion  with  respect  to  hydrogen 

=  ionic  charge 

uncertainty  factor  explained  in  section  5 

l/X^l/3,  fraction  of  space  occupied  by  matter,  explained  under 
Results  (5). 

Evaluating  equation  30  for  the  case  of  a  constant  temperature  of 
6  o 

T  =  10°  K  we  arrive  at  the  following  expression: 


(30) 
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tot 


v-23 


N  ^UX 
e 


2.  25  X  10  erg  cm  /sec 


(31) 


With  the  various  ionic  abundance  as  calculated  for  T  =  1  x  10^  °K  but 
with  T  =  7.  5  X  10^  °K  we  obtain 


tot 


N  ^LTX 
e 


=  2.  34  X  10 


■23  o„  3/ 

K  erg  cm  Aec 


(32) 


The  individual  sources  of  radiative  loss  have  been  computed 
separately,  and  the  resiilts  are  shown  in  Table  3.  The  total  has  been  checked 
by  substituting  equation  (32)  into  equation  (1). 


4.  Results 

In  summarizing  the  radiative  losses  due  to  the  various  processes 
described  above  (see  Table  3)  we  make  the  following  observations: 

(1)  The  overwhelming  contribution  to  the  radiative  loss  occurs  within  15 
solar  radii.  As  a  result  the  energy  losses  differ  very  slightly  when 
different  models  of  temperature  and  electron  density  are  employed,  since 
the  models  differ  significantly  from  one  another  only  at  large  distances 
from  the  sun. 

(2)  The  single  most  important  source  is  the  X-ray  lines  which  contributes 
a  little  over  half  the  total. 

(3)  The  heavy  ion  which  makes  the  largest  contribution  to  the  radiative 
loss  in  the  ultraviolet  region  is  Si  XI.  However,  it  is  probably  not  safe 
to  give  strong  weight  to  the  calculations  for  individual  lines  or  ions. 

(4)  The  total  radiative  loss  from  the  heavy  elements  in  the  ultraviolet  region 

21 

is  about  9x10  ergs/sec.  However,  in  the  calculations  only  the  transitions 
from  the  first  excited  level  to  the  ground  level  were  accovmted  for. 
Transitions  between  the  higher  levels,  with  subsequent  cascades  to  the 
gro\md  level,  were  not  calculated  since  no  reliable  f-values  are  available 
for  such  transitions.  However,  in  order  to  estimate  the  relS’tive  importance 
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of  the  higher  levels,  we  may  look  back  at  the  calculations  for  H  and  He. 

Here  the  calculations  did  take  into  account  the  upper-level  transitions  and 

subsequent  cascade  to  the  ground  level;  such  transitions  constituted  about 

1/3  of  the  bo\ind-  bound  radiative  loss  for  H  I  and  He  II.  If  it  is  also  true  that 

1/3  of  the  radiative  loss  from  bound-bound  transitions  in  heavy  elements  comes 

from  the  upper  levels,  then  the  total  radiative  loss  in  the  ultraviolet  from  the 

21  21 

heavy  elements  should  be  13.  5  x  10  ergs/sec  instead  of  9  x  10  ergs  sec. 

22 

(5)  The  total,  computed  radiative  loss  from  all  sources  is  ^v-’4.  9  x  10 

2  2 

ergs/sec.  However  we  must  note  here  that  actually  <  Ne  >  =  1/X'^Ne  ^ 

where  1/X  can  be  thought  of  as  the  fraction  of  space  occupied  by  matter 

(see  Unsold  I960  and  Brandt  1962b).  Thus  we  are  underestimating  the 
2 

quantity  ,  which  appears  in  every  expression  which  we  evaluate  by  the 

factor  1/X.  The  estimates  of  the  quantity  1/X  (valid  near  5R  )  are  about 

1/5  (Unsold  I960).  This  correction  would  raise  the  radiative  loss  to  about 
23 

2  X  10  ergs/sec. 

(6)  The  diffuse  intensity  of  the  permitted  ultraviolet  lines  contributing  to 
the  radiative  energy  loss  is  very  small  (47rilc  R)  for  observations  at 
elongations  corresponding  to  distance  from  the  center  of  the  sun  of  4R^  or 
greater. 

5.  Comparison  with  Other  Data 


To  check  the  calculations  we  may  compare  them  with  the  rocket 

2 

observations  of  the  solar  flux  per  cm  in  the  ultraviolet  and  X-ray 
regions  at  the  earth's  distance  from  the  sun.  To  carry  out  the  comparison 
we  integrate  eqiiation  (1)  backwards  to  the  inner  boundary  of  the  corona 
(r  =  1.  015  R  )  using  the  electron  densities  for  solar  minimum  given  by 
de  Jcvger  (1959).  This  is  essentially  a  calculation  of  J  d V ;  we  note 

that  the  same  integral  for  the  solar  maximum  data  given  by  de  Jager 
is  16  times  the  value  for  solar  minimum.  The  comparison  can  at  best  be 
a  rough  one,  principally  because  (1)  the  experimental  results  are 
tincertain  by  a  factor  of  3  to  5  and  (2)  and  because  of  the  uncertain  role  of 
the  transition  region  in  the  production  of  vdtraviolet  photons. 


Elwert  (1961)  quotes  the  measured  values  0„  13,  0.  3,  and  1.  0 

2 

erg/ cm  sec  for  the  minimum,  mean,  and  maximum  solar  flux  outside  the 
earth's  atmosphere  in  the  X-ray  region  between  44  and  100  His 
corresponding  theoretical  calculations  give  a  mean  value  of  0.12  XJ  erg/cm' 
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sec,  wheretJ  indicates  the  uncertainties  in  the  theoretical  cross  sections. 

Elwert's  calculation  agrees  with  the  observations  if  U  is  taken  between 

2  and  4.  Our  calculations,  on  the  other  hand,  give  a  result  of  0.  05  "O' 

2 

erg/cm  sec;  for  solar  minimum.  We  consider  that  only  one-half  of  the 
X‘<ray8  produced  escape  into  space,  and  that  X^2.  Our  calculations  then 
require  2.  Such  a  value  is  quite  reasonable  since  U  must  include 

errors  in  the  observational  results  as  well  as  in  the  theoretical  cross 
sections,  it  also  must  include  all  sources  of  radiative  loss  which  were  not 
taken  into  account  in  the  calculations  (incompleteness). 

Hinterreger  (1961)  has  also  published  observations  covering  the 

ultraviolet  and  X-ray  regions.  His  total  observed  flux  due  to  the 

metallic  lines  in  the  region  X  <  1300  A. ,  corresponding  in  the  most  part 

2 

to  our  calculated  lines,  is  about  1.  5  erg/cm  sec.  Our  theoretical  result 

2 

(assuming  X^2  and  that  half  the  photons  escape)  is  0.  OSUerg/cm  sec 

2 

for  solar  minimum  and  1.  SCJerg/cm  -sec  for  solar  maximum.  Since 
Hinterreger's  observations  refer  to  solar  maximvun  conditions  (Cragg  1961), 
we  haveU”^ls  1,  i.  e.  ,  satisfactory  agreement. 

Woolley  and  Allen,  in  their  1948  estimate  of  radiative  losses  from 
the  corona,  obtain  a  value  for  line  emission  in  the  ultraviolet  which  is 
about  10  times  the  value  we  have  obtained  for  solar  minimum.  This  figure 
is  easily  explained  by  the  following  two  considerations:  (1)  use  of  the 
mean  Baumbach- Allen  densities  raises  their  value  by  a  factor  of  three 
and  (2)  lacking  both  observations  and  f-values  for  the  far  ultraviolet,  they 
assumed  all  line  emissions  to  be  concentrated  into  ultimate  lines  for  which 
f  =  1,  whereas  the  average  f-value  for  lines  in  this  region  is  much  less  than 
this  (Varsavsky  1961).  Our  calculations  are  also  in  satisfactory  agreement 
with  the  recent  calculations  of  Allen  (I960),  when  differences  in  assumed 
abundances  and  f-values  are  considered.  We  expect  that  our  calculations 
are  an  improvement  over  Allen's  (1961)  because  of  the  use  of  Varsavsky's 
(1961)  f  -values. 
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6.  Consequences 

Brandt  and  Michie  (1962,  eq.  20)  derive  an  expression  relating 
the  conduction,  velocity  space  anisotropy,  expansion  velocity,  temperature, 
and  heliocentric  distance  at  points  in  the  corona.  The  derivation  of  this 
expression  is  based  on  the  assumption  that  ^  e/’^m'  ratio  of 
energy  flow  to  mass  flow  across  a  spherical  surface  at  distance  r  from 
the  sun,  is  a  constant,  independent  of  r,  beyond  a  certain  distance  which 
is  assumed  to  be  4  solar  radii.  The  derivation  did  not  take  into  account 
radiative  losses.  We  now  re-derive  the  expression  with  the  inclusion  of 
a  radiative  loss  term  in  order  to  demonstrate  that  the  assumption  made 
earlier  still  is,  in  fact,  valid. 

Equation  (13)  of  Brandt  an 

m 


modified  to  become 


e  * 


d  Michie  (1962)  for  the  energy  flow  is 


2iTv^dpdv  +  2Nw<v 


(33) 


where  F(?,  (i  )  is  the  distribution  function  in  terms  of  the  peculiar 
velocity,  v  is  the  peculiar  velocity  vector,  w  is  the  mean  expansion 
velocity,  6  is  the  angle  between  the  radius  vector  and  the  peculiar 
velocity,  tr,  p  =  cos  ^  ,  m  =|l/^mpj  (the  effective  particle  mass). 


N=  J  •  2uv^dpdv  ,  (34) 


and  is  tke  radiative  flux  (The  subscript  r  refers  to  the 

radial  component.  )  Since  the  mass  flow  is  given  by 
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m 


4i!r 


mNw 


(35) 


we  have 


1 


jiv^F 


Zw^djjwlv  +  2  < 


(36) 


+  w 


2 


^rad 

+  - 

mNw 


Let 


3kT 


m 


(37) 


(38) 


(39) 


and  (36)  modifies  to 


^  e 


2GM 

_ • 

r 


m  )  mNw 


(40) 
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Evaluating  the  quantity  at  some  inner  position  we  obtain 


I  w 
o  ^ _ o_ 


GM  kT  V 

_ ®  +  _ 2T2.—  + 


m 


^  rad, 


mN  w 
o  o 


GM 


kTy  1 

+  -  + 


m 


mNw 


rad,  r 


(41) 


below  r  and  d? 


where  ,  is  the  radiative  flux  originating  -  , 

rad,  or  ®  °  o  rad,  Tq  r 

is  the  radiative  flux  originating  between  r^  and  r.  Noting  that 


rad,  o 


mN  w 
o  o 


rad,  or 
mNw 


and  simplifying  the  resulting  expression,  we  obtain 


(42) 


kTy  7  j  kT  jr 

'  =  (l/2)(I^-I+w^^-w'‘)  -  GM^(l/r^-l/r)  +  ° 


ir 


m 


m 


4TimC  l^r 


e  (r)4irr  dr 


where  we  have  written 


(43) 


C^rad, 


lir  t/*  « (r)r^ 


dr 


rpr 


mNw 


4irr  Nmw 


4nmC 


^  • 


4firr  dr  (44) 


In  an  entirely  analogous  manner  we  may  include  a  term  which  accounts  for 
the  energy  deposition,  C.  (r),  between  r^  and  r,  giving  as  a  final  expression 
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2  2 

— 1-  =  (l/2){I^-I+w/-w^) 


kT  y 

GM  (l/r  'l/r)  + - - — 2 —  + 

«  o 


m 


4iTr^dr 


(45) 


Evaluation  of  the  quantities  in  eqviation  (45),  taking  T  =  10^  °K  and  ^  =  2.  5, 


shows  that 
of 


kTY 

m 


13 

5  X  10  while  the  radiative  loss  term  has  a  value 


12  72 

2  X  10  if  the  flux  near  the  earth  is  taken  as  3  x  10  electrons/cm  -sec. 


The  quantity  y,  which  indicates  the  velocity  space  anisotropy,  must  satisfy 
the  relationship  (Brandt  and  Michie  1962  eq.  (21))  3/2  ^  “X  ^  9/2.  Therefore 
for  all  temperatures  greater  than  10^  °K  we  can  say  that  the  radiative. loss 
term  is  negligibly  small  in  comparison  with  the  quantity  .  Even 


m 


if  the  radiative  loss  must  be  corrected  by  a  factor  3  this  statement  still 

At  o 

holds.  Only  when  the  temperature  drops  to  10  K  or  less  does  the 
radiative  loss  amount  to  an  appreciable  fraction  of  ;  however 


temperatures  as  low  as  this  do  not  occur  in  the  models  of  the  interplanetary 

5  o 

gas  which  we  use,  and,  in  fact,  a  temperature  of  about  2x10  K  was 
found  from  the  Mariner-II  results  by  Neugebauer  and  Snyder  (1962).  Thus 
we  may  safely  assume  that  is  a  constant  beyond  r  =  4R^;  this 

conclusion  may  be  regarded  as  the  most  important  consequence  of  the 
radiative  loss  calculation.  It  is  still  subject  to  the  condition,  ho^^evej:,  that 
the  energy  deposition  also  be  small  beyond  4R^.  Further  calcvilation  is 
necessary  in  order  to  bear  out  this  hypothesis,  but  if  A  (r)  is  small  it 
serves  to  cancel  the  effects  of  c  (r)  and  therefore  to  reinforce  the  conclusion. 


The  calculation  also  has  implications  regarding  the  quantity  I^. 

The  expression  mw^NI^  represents  the  energy  flux  transported  by  conduction, 

evaluated  at  the  point  r  =  4R  .  This  was  calculated  by  Brandt  and  Michie; 

*  2  2 

they  obtained  a  value  of  4  x  10  ergs/^m  which  they  took  to  be  a  maximum 
value,  since  it  included  a  number  of  quantities  which  had  not  been 
specifically  taken  into  account,  notably,  the  radiative  losses.  To  obtain 
consistency  of  the  theoretical  model  with  observation,  they  found  that  had 
to  be  reduced.  This  reduction  can  be  accounted  for  by  one  or  more  of  the 
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following  effects:  (1)  coronal  magnetic  fields  which  would  tend  to  inhibit 

conduction:  (2)  radiative  losses  beyond  4R  ;  and  (3)  an  actual  particle  flux 

in  the  solar  wind  greater  than  that  adopted  on  the  empirical  model.  As  a 

result  of  the  present  calculation,  the  hypothesis  that  radiative  losses 

between  4R^  and  the  orbit  of  the  earth  form  a  significant  part  of  must  be 

discarded.  This  result  is  in  agreement  with  the  recent  rocket  results 

8  Z 

(Neugebauer  and  Snyder  1962)  which  give  a  flux  ^1  x  10  electrons/cm  -sec. 
Thus,  all  or  part  of  the  discrepancy  is  resolved  with  the  new  value  for  the 
flux^point  (3)  above^.  We  still  have  the  possible  problem  of  the  coronal 
magnetic  field. 

The  existence  of  irregular  or  non- radial  magnetic  fields  could 
easily  be  quite  important  in  determining  the  thermal  conductivity  in  the 

corona  near  4R  .  The  thermal  conductivity  appears  to  be  reduced  by  a 

®  2  2-1 
transverse  magnetic  field  by  a  factor  of  (1  +  t^^  )  (Spitzer  1956, 

Cowling  1954),  where  is  the  cyclotron  frequency  for  electrons  and 

t^^  is  the  collision  time  for  electrons.  Thus,  for  the  magnetic  field  not 

to  be  important  in  inhibiting  heat  conduction,  we  must  have 


eB 


cm 

e 


^3/2 

-  <  1 

15N 

e 


(46) 


When  values  of  the  various  parameters  appropriate  to  the  corona  at  4R 

W 

are  inserted  into  equation  (46),  we  find  that  a  transverse  magnetic  field 
must  be  less  than  10  gauss  not  to  inhibit  conduction.  Since  the 
magnetic  field  expected  in  the  neighborhood  of  4R  is  many  orders  of 
magnitude  greater  than  this  limit,  we  conclude  that  any  reasonable 
transverse  magnetic  field  at  4R^  will  seriously  inhibit  conduction.  Thus, 
a  knowledge  of  the  structure  of  the  solar  magnetic  field  at  some  distance 
from  the  s\in  would  seem  to  be  most  desirable.  We  note  that  the 
contribution  to  the  radiative  loss  teitm  from  magnetic  bremmstrahlung 
is  probably  quite  small. 

It  should  be  mentioned  that  the  low  densities  found  on  the  semi- 
empirical  model  (Brandt  and  Michie  1962)  and  the  Mariner- II  results 
(Neugebauer  and  Snyder  1962)  are  supported  by  additional  evidence.  The 
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low  densities  are  favored  on  the  basis  of  Wurm's  (1962)  study  of  comet 

«l 

comas  and  on  the  basis  of  Opik's  (1956)  estimate  of  the  upper  limit 

3  ‘ 

(  ^  5  electrons/cm  )  of  the  density  of  the  solar  wind  consistent  with  the 
existence  of  meteoric  particles  in  the  interplanetary  medium.  The  plasma 
probe  measures  of  Bonetti,  Bridge,  Lazarus,  Lyon,  Rossi  and  Scherb 
(1962)  are  also  consistent  with  a  low  mean  density  (see  Piddington  1962). 

Finally,  we  remark  that  equation  (30)  is  a  general  expression 
which  can  be  used  to  evaluate  the  radiative  losses  in  any  rarefied  gas  with 
a  temperature  on  the  order  of  10^  °K  in  which  statistical  equilibrium  and 
the  Maxwellian  velocity  distribution  hold.  In  the  paper  the  expression  was 
evaluated  over  the  ranges  of  temperature  and  pressure  in  the  outer  solar 
corona;  however  it  could  also  be  applied  in  other  situations.  In  particular, 
the  spectra  of  the  recurrent  nova  RS  Ophiuchi  showed  coronal  forbidden 
emission  lines  in  the  outbursts  of  1933  and  1958,  indicating  a  temperature 
^  10^  °K.  Wallerstein  (1961)  has  estimated  the  electron  density  after  the 
outburst  in  the  circumstellar  region  to  be  on  the  same  order  as  that  of  the 
inner  corona  of  the  sun,  indicating  that  the  work  presented  here  could  also 
be  of  value  in  the  explanation  of  the  coronal  lines  of  RS  Ophiuchi. 
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TABLE  1 


Abundances  of  Various  Elements 


Element 

log  Ne/Nh  +  12 

Element 

log  Ng/Njj  +  12 

H 

12.0 

Si 

7.  5 

He 

11.0 

S 

7.  3 

C 

8.7 

Ca 

6.  2 

N 

8.  0 

Fe 

6.  6 

O 

9.0 

Ni 

5.9 

Na 

6.3 

Ne 

8.  1 

Mg 

7.4 

A 

6.  3 

A1 

6.  2 

TABLE  2 


Ion  Densities  and  Bovind-bound  Radiative  Losses  for  Heavy  Elements 


Ion 

Pi  =  n./Nh 

>i(A) 

Radiative  Loss  x 
10”^^  erg/sec. 

X-Ray  Lines  Included 

in  Calculations 

Mg  IX 

1.  5  X  10 

368.  1 

.  25 

Element 

All) 

Mg  X 

7.5  X  10  5 

625.  3,  609.  8 

.  55  — 

Mg  XI 

1..6x  lO'J 

-  -  -  + 

.  .  - 

C  IV 

312 

O  VI 

1.4  X  10"? 

1037.6,  1031.9 

.  24 

N  V 

210 

O  VII 

9.7  X  10'!: 

21.  60 

•  •  . 

N  VI 

29 

Ni  XII 

1. 1  X  10"'_ 

299,  322 

.  06 

N  VII 

19 

Ni  XIII 

2.86  X  10, 

304-332 

.  10 

O  VI 

150 

Ni  XIV 

2.7  X  10‘, 

311-323 

.  05 

O  VII 

22 

Ni  XV 

1.0  X  10‘, 

316,  333,  346 

.02 

Ne  VII 

95 

A  DC 

8.0  X  10', 

... 

Ne  Vm 

85 

A  X 

9.0  X  10, 

166,  171 

.  27 

Mg  VIII 

75 

A  XI 

3.  2  X  10"^ 

190 

.  18 

Mg  DC 

63 

C  V 

3.0  X  10^ 

40.  27 

.  21 

Mg  X 

58 

C  VI 

2.5  X  10"^ 

33.  2 

.26 

Si  vni 

67 

Na  IX 

2  X  10*  7  , 

681.  7,  694.  3 

.02 

Si  IX 

55 

Na  X 

1.8  X  10'^ 

- + 

... 

Si  X 

51 

Ne  VIII 

2. 5  X  10“4 

792.  5,  782,  4 

.  25 

Si  XI 

42 

Ne  IX 

1.2  X  10 

---  + 

... 

S  VII 

70 

A1  X 

4.  8  X  10*, 

332.  9 

.08 

S  VIII 

60 

A1  XI 

7.0  X  10*5 

550.  0,  568.  5 

.  04 

S  DC 

52 

Si  IX 

4.7  X  lo'r 

342- 350 

.42 

S  X 

43 

Si  VIII 

4.7  X  10*'. 

314,  316,  320 

.  04 

Fe  DC 

150 

Si  X 

1.38  X  10'^ 

356.  1,  347.4 

1.  05 

Fe  X 

97 

Si  XI 

1.38  X  10*^ 

303.  6 

1.  76 

Fe  XI 

88 

Si  XII 

4.  5  X  io*r 

521.  1,  499.  3 

.  27 

Fe  XII 

73 

S  VIII 

6.  3  X  10  , 

198.  6,  202.  6 

.  14 

Fe  XIII 

65 

S  IX 

4.0  X  10"? 

222-228 

.56 

Fe  XIV 

62 

S  X 

9.8  X  10*2 

258,  260,  264 

.64 

Fe  XV 

55 

S  XI 

6.2  X  io*r 

290 

.  36 

S  XII 

6.  0  X  10  , 

302,  290 

.03 

Fe  XII 

3.  5  X  10*, 

356,  360,  369 

.  07 

Fe  XIII 

9. 6  X  10"5 

360,  373,  386 

.20 

Fe  XIV 

1. 8  X  10"? 

346,  370 

.  37 

Fe  XV 

1.2  X  10', 

288 

.44 

Fe  XVI 

2.6  X  10"^ 

361.7,  336.2 

.04 

N  VI 

7.5  X  10"? 

- + 

... 

Cai  XI 

1.  5  X  10*° 

-  -  - 

TOTAL 

9. 0  X 10^^  erg/sec 

26.  2  X  10^^ 

erg/sec 

no  lines  calculated  for  this  element 
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TABLE  3 

Contributions  of  the  Various  Processes  to  the  Radiative  Loss  (equation  1) 

(Unit:  10^^  erg/sec) 


Hydrogen,  total 

3.  5 

Free-free 

1.8 

Free-bound 

1.  1 

Bound-bound 

0.  6 

(upper  limit) 

Helium,  total 

4.  3 

Free-free 

0.7 

Free-bound 

1.  8 

B  ound-  bo\md 

1.  8 

(upper  limit) 

Heavy  elements,  total 

41.  8 

Free-free 

0.  1 

Free-bo\md 

1.9 

Forbidden  lines 

0.  1 

Permitted  lines+ 

9.  0 

(upper  levels)* 

4.  5 

X-ray  lines 

26.2 

Total 

49.  6 

+For  details  see  Table  2 
*See  (4)  under  results. 
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